library(dplyr)
library(ggplot2)
library(plotly)
theme_set(theme_bw())
# data ----
source(here::here('r', 'bridge_data.r'))
data.frame(age = ages,
year = readRDS(here::here('data', 'rtmb', 'report1.rds'))$fish_age_pred,
groups = 'rtmb') -> df
names(df) <- c('age', fish_age_yrs, 'groups')
read.csv(here::here('data', 'admb', 'fac.csv')) %>%
mutate(groups = ifelse(groups=='pred', 'admb', groups)) %>%
bind_rows(df %>%
tidyr::pivot_longer(-c(age, groups), names_to = 'year') %>%
mutate(year = as.numeric(year))) %>%
ggplot(aes(age, value, color = groups, linetype=groups)) +
geom_point() +
geom_line() +
facet_wrap(~year) +
scico::scale_color_scico_d(palette = 'roma', end = 0.8) +
scale_linetype_manual(values = c(0,1,1), guide = guide_none()) -> f1
ggplotly(f1)